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ABSTRACT 

Prestellar cores are unique laboratories for studies of the chemical and physical conditions 
preceding star formation. We observed the prestellar core L1544 in the fundamental transition of 
ortho-H 2 D + (li,o-li,i) at different positions over 100", and found a strong correlation between its 
abundance and the CO depletion factor. We also present a tentative detection of the fundamental 
transition of para-D 2 H + (li o — lo,i) a * the dust emission peak. Maps in N 2 H + , N 2 D + , HCO + 
and DCO + are used, and interpreted with the aid of a spherically symmetric chemical model 
that predicts the column densities and abundances of these species as a function of radius. The 
correlation between the observed deuterium fractionation of H3" 1 ", N 2 H + and HCO + and the 
observed integrated CO depletion factor across the core can be reproduced by this chemical 
model. In addition a simpler model is used to study the H 2 D + ortho-to-para ratio. We conclude 
that, in order to reproduce the observed ortho-H 2 D + observations, the grain radius should be 
larger than 0.3 /im. 

Subject headings: ISM: molecules — ISM: individual (L1544) — radio lines: ISM 



1. Introduction 

Deuterium bearing species are good probes of 
the cold phases of molecular clouds prior to star 
formation. Many recent observations point to 
the fact that their abundances relative to their 
fully hydrogcnatcd forms are larger, by factors up 
to 10 5 , than the solar neighborhood value of ~ 
1.5 x 1CT 5 found by Linsky (2003). The deu- 
terium fractionation has been evaluated in prestel- 
lar cores and low mass protostars from observa- 
tions of HCO+ and N 2 H+ (Butncr et al. 1995; 
Williams et al. 1998; Caselli et al. 2002a; Crapsi 
ct al. 2004, Crapsi et al. 2005), H 2 CO (Loinard 
et al. 2001; Bacmann et al. 2003), H 2 S (Vastel 
ct al. 2003), HNC (Hirota et al. 2003), CH 3 OH 
(Parise et al. 2004), and NH 3 (Roucff ct al. 2000, 
Tine et al. 2000). The chemical fractionation 
process in the gas-phase mainly arises from the 
difference between the zero-point energies of H 2 
and HD. Almost incredibly, this can lead to a de- 
tectable quantity of triply deuterated molecules 
like ND 3 (Lis et al. 2002; van der Tak ct al. 
2002) and CD 3 OH (Parise et al. 2004). The lat- 
ter is thought to be formed mainly on dust grain 
surfaces (Charnley et al. 1997) in regions where 
the gas-phase [D]/[H] ratio is enhanced to values 
larger than ^0.1 (Caselli et al. 2002c), as in the 
cold cores (see below). In molecular clouds, hy- 
drogen and deuterium are predominantly in the 
form of H 2 and HD respectively. So the HD/H 2 
ratio should closely equal the D/H ratio. Since the 
zero-point energies of HD and H 2 differ by ~ 410 
K, the chemical fractionation will favor the pro- 
duction of HD compared to H 2 . In the dense, cold 
regions of the interstellar medium (T ~ 10 K), D 
will be initially nearly all absorbed into HD. The 
abundant ion available for interaction is Hg , which 
gives H 2 D+: 

H++HD< — ► H 2 D+ + H 2 + 230K (1) 

The reverse reaction does not occur efficiently in 
the cold dense clouds where low-mass stars form, 
and where the kinetic temperature is always be- 
low 25 K, the "critical" temperature above which 
reaction (1) starts to proceed from right to left 
and limits the deuteration. Therefore, the degree 
of fractionation of H 2 D + becomes non- negligible. 
This primary fractionation can then give rise to 
other fractionations and form D 2 H+ and D^ as 



first suggested by Phillips & Vastel (2003): 

H 2 D+ + HD < — > D 2 H+ + H 2 + 180K (2) 

D 2 H+ + HD < — ► D+ + H 2 + 230K (3) 

We present in Figure 1 the main reactions involv- 
ing these molecules. Note that the effect of the 
recombination of H^ with electrons in the gas is 
negligible because of the low electron density in 
such regions. However, the effect of recombination 
with electrons on negatively charged grain surfaces 
becomes important when depletion increases (cf. 
Walmsley et al. 2004). The dissociation of the 
deuterated forms of H^ is then responsible for the 
enhancement in the [D]/[H] ratio. One specific 
parameter can enhance this process: the deple- 
tion of neutral species (in particular, the abun- 
dant CO) from the gas-phase (cf. Dalgarno & 
Lepp 1984). In fact, the removal of species which 
would normally destroy H^ (e.g. CO; Roberts 
& Millar 2000) means that H^ is more likely to 
react with HD and produce H 2 D + , D 2 H + and 
Dg". The first model including D 2 H + and D^ 
(Roberts, Herbst & Millar 2003) predicted that 
these molecules should be as abundant as H 2 D+ 
(see also Flower et al. 2004). 

Gas phase species are expected to be depleted 
at the center of cold, dark clouds, since they tend 
to stick onto the dust grains. A series of recent 
observations has shown that, in some cases, the 
abundance of molecules like CO decreases toward 
the core center of cold (< 10 K), dense (> a few 
10 4 cm" 3 ) clouds. L1498: Willacy et al. (1998), 
Tafalla et al. (2002, 2004); IC 5146: Kramer et al. 
(1999), Bergin et al. (2001); L977: Alves et al. 
(1999); L1544: Caselli ct al. (1999), Tafalla et 
al. (2002); L1689B: Jessop & Ward-Thompson 

(2001) ; Bacmann et al. (2002); Redman et al. 

(2002) ; B68: Bergin et al. (2002); L1517B: Tafalla 
et al. (2002, 2004); L1512: Lee et al. (2003); Oph 
D: Bacmann et al. (2003), Crapsi et al. (2005); 
L1521F: Crapsi et al. (2004); L183 (L134N): Pa- 
gani et al. (2005). These decreases in abundance 
have been interpreted as resulting from the deple- 
tion of molecules onto dust grains (see, e.g., Bergin 
& Langer 1997; Charnley 1997). It is now clear 
that these drops in abundance are typical of the 
majority of dense cores (see Tafalla & Santiago 
2004 for the case of L1521E, the first starless core 
to be found with no molecular depletion). 
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In one of the most heavily CO depleted prestel- 
lar cores, L1544, Casein et al. (2003) detected 
a strong (brightness temperature ~ 1 K) ortho- 
H2D + (loi-ln) line, and concluded that H 2 D + is 
one of the main molecular ions in the central re- 
gion of this core. Encouraged by laboratory mea- 
surements (Hirao & Amano 2003), Vastel et al. 
(2004) detected the para-D 2 H+ molecule in its 
ground transition at 692 GHz. They found that, 
in the prestellar core 16293E, D2H + is as abun- 
dant as H 2 D + . These studies supported chemical 
modelling and the inclusion of multiply deuterated 
species (Roberts et al. 2003; Walmsley et al. 2004; 
Roberts et al. 2004; Flower et al. 2005; Aikawa et 
al. 2005). It appears that in dark clouds affected 
by molecular depletion, the deuterated forms of 
the molecular ion Hjj~ are unique tracers of the 
core nucleus, the future stellar cradle. Thus, their 
study becomes fundamental to the unveiling of 
the initial conditions of the process of star forma- 
tion (kinematics, physical and chemical structure 
of pre-stellar cores). 

In this paper, we present new observations of 
the H 2 D+ (lio-ln) line toward L1544, mapped 
over ~ 100" of the dust peak emission, as well as 
a tentative detection of the D 2 H+ (lio-loi) line. 
Caselli et al. (2003) roughly estimated the size of 
the H 2 D + emitting region in this pre-stellar core 
and suggested a radius of about 3000 AU. But 
this was based on a five-point map and cannot 
put stringent constraints on the chemical struc- 
ture. Also a parallel study has been done by van 
der Tak et al. (2005) with the analysis of the line 
shape profile. Here, the H 2 D + map is also com- 
pared with other high density tracer maps. Due 
to the poor atmospheric transmission at the fre- 
quency of the para-D 2 H + fundamental line, this 
study is limited to the ortho-H 2 D + fundamental 
line. We will present, in the last section, the per- 
spectives on this work that will be opened up by 
future observatories. 

2. Observations 

The observations were carried out at the Cal- 
tech Submillimctcr Observatory (CSO), between 
November 2003 and February 2005, under good 
weather conditions (225 GHz zenith opacity al- 
ways less than 0.06) where the atmospheric 
transmission is about 30% at 372 GHz and less 



than 20% at 692 GHz. Scans were taken to- 
ward the peak of the 1.3 mm continuum dust 
emission of L1544 (a 200 o = 05 h 04 m 17 s . 21, 
^2000 = + 25° 10' 42.8") using the chopping 
secondary with a throw of 3'. The 345 GHz (re- 
spectively 650 GHz) sidecab receiver with a 50 
MHz acousto-optical spectrometer backend was 
used for all observations with a velocity resolu- 
tion of 0.06 km s _1 (respectively 0.03 km s _1 ) 
i.e. ~ 1.6 channels. At the observed frequen- 
cies of 372.421385(10) GHz for the H 2 D+ (1 10 - 
ln) and 691.660483(20) for the D 2 H+ (li -l i) 
lines (Amano & Hirao 2005), the CSO 10.4 meters 
antenna has a HPBW of about 20" and 11" re- 
spectively. We mapped the area in H 2 D+ around 
the dust peak position with a grid spacing of 20" 
and used the value at the peak from Caselli et 
al. (2003) and integrated longer. The beam ef- 
ficiency at 372 GHz (respectively 692 GHz) was 
measured on Venus, Saturn and Jupiter and found 
to be ~ 60% (respectively ~ 40%). Pointing was 
monitored every 1 hour and half and found to be 
better than 3". In the case of H 2 D + , the emission 
is extended compared to the beam size of CSO: 
the efficiency is then about 70%. If the emission 
in D 2 H + is extended compared to the beamsizc 
of 11", then the efficiency at 692 GHz is about 
60%. The data reduction was performed using 
the CLASS program of the GAG software devel- 
oped at IRAM and the Observatoire de Grenoble. 

Figure 2 shows the H 2 D+ and D 2 H + spectra 
observed toward the dust peak position. A Gaus- 
sian fit to the H 2 D + line gives a LSR velocity of 
^7.3 km s _1 , but two peaks are clearly visible, 
with velocities 7.1 and 7.4 km s _1 (van der Tak 
et al. 2005), as also seen in other tracers (Tafalla 
et al. 1998; Caselli et al. 2002a). This central 
position has been originally observed by Caselli et 
al. (2003) and studied in detail by van der Tak 
et al. (2005); here, we improved the sensitivity 
and used the new value of the ortho-H 2 D + (li - 
ln) line frequency, recently measured by Amano 
& Hirao (2005). The on-source integration time 
for D 2 H+ is about 230 min. The D 2 H+ feature 
can be fitted with a gaussian with T a *=0.30 ± 
0.07 K, Av=0.08 ± 0.04km s" 1 and V LS fl=7.29 ± 
0.03 km s _1 . The solid vertical line corresponds to 
the velocity from this gaussian fit. It is consistent 
with the central velocity of the H 2 D + line of 7.28 
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± 0.06 km s _1 . However, it is difficult to believe 
that the observed linewidth represents the real 
linewidth of the transition as the expected ther- 
mal linewidth is about 0.25 km s _1 for a kinetic 
temperature of 7 K, as predicted by dust temper- 
ature measurement (Evans et al. 2001; Zucconi et 
al. 2001). This is about 3 times larger than the 
observed linewidth. Formaldehyde observations of 
Bacmann et al. (2002) suggest larger gas temper- 
atures (up to 9 K), but H2CO is likely frozen in 
the core center, so that the measured tempera- 
ture probably reflects the warmer core envelope 
(Young et al. 2004). At high densities (higher 
than - 3 10 4 cm~ 3 ), Young et al. (2004) find the 
gas and dust temperature to be between ~ 7 and 
9 K, consistent with what was found by Tafalla 
et al. (2002) using ammonia. We will use in the 
following a temperature of 8 K in this cloud. The 
signal to noise ratio of our D 2 H + observations is 
not sufficient to get constraints on the kinemat- 
ics of the source and the fitted linewidth probably 
does not represent the real linewidth. In absence 
of a possible explanation for the narrow linewidth 
for the D 2 H+ profile, we will consider, in the fol- 
lowing, that we only have a tentative detection, 
and will then use an upper limit. 

Figure 3 shows the H 2 D+ spectra around the 
central position (0",0") of the dust peak emission. 
The offset positions shown in the upper left are in 
arcseconds. A Gaussian fit for each detected line 
is plotted, using the CLASS program, and the line 
parameters are presented in Table 1. A double- 
peak profile seems to appear in the south-east as 
well as in the central part of the map. Consider- 
ing the rms value, it is only possible to say that 
this non Gaussian profile is localized in the cen- 
tral positions around the dust peak emission. A 
possible interpretation could be the presence of 
two different layers with different velocities along 
the line of sight. Another explanation could be 
that the observed profiles are affected by absorp- 
tion in a low-density (10 4 cm -3 ) foreground layer 
redshifted (~ 0.08 km s _1 ) relative to the high- 
density core, as found by (Williams et al. 1999) 
in the case of N 2 H + (l-0) mapped at high spatial 
resolution. A more detailed study of the H 2 D + 
line profile toward the L1544 dust peak has been 
recently carried out by van der Tak et al. (2005), 
who suggested that the observed H 2 D + line, be- 
sides being broadened by the central infall, can 



also be absorbed in the outer parts of the core. 
The presence of a central dip in the H 2 D+ profile 
of at least four spectra across the L1544 map (see 
Fig. 3) favours this scenario. 

Figure 4 shows the integrated intensity map 
(/ T m bdv) of ortho-H 2 D+ (I10-I11), together with 
maps of N 2 H+ (1-0) and N 2 D+ (2-1) obtained by 
Casein et al. (2002a) and the 1.3 mm continuum 
emission map from Ward-Thompson et al. (1999), 
smoothed to a resolution of 22". We note the close 
similarity between the H 2 D+ and the N 2 D+ maps 
and this will be discussed in the next sections. 
In this paper, we studied the chemistry using the 
maps made by Caselli et al. (2002a) in H 13 CO+, 
HC 18 0+, DCO+, D 13 CO+, C 17 and C 18 with 
the IRAM 30 m telescope. 

3. Column density and abundance deter- 
minations 

The observed molecular ions maps presented in 
Figure 4 show a general correlation, despite dif- 
ferent beamwidths, with the distribution of dust 
continuum emission, in contrast to C 18 (1-0) and 
C 17 (1-0) (Caselli et al. 1999), which give clear 
evidence for depletion of CO at positions close to 
the continuum peak. H 2 D+ (lro-lii), N 2 D+ (2-1) 
and to a lesser spatial extent N 2 H + (1-0) appear 
to trace the dust continuum. From these maps 
N 2 n + does not seem to be depleted at the dust 
peak position. 

In order to compare the observed species and put 
constraints on chemical models, we need to in- 
fer the column densities and abundances of H 2 D + 
and D 2 H+ (defined as A^(i)/A(H 2 ) for a generic 
species i, with A(H 2 ) derived from the 1.3 mm 
dust continuum emission map of Ward-Thompson 
et al. 1999). Assuming LTE conditions, we can es- 
timate the optical depth at the line center from the 
observed line intensities: 

T mb = [J,(T ex )-J ! ,(T bg )](l-e- T ) (4) 

where J V {T) = {hv/k)/{e hv / kT - 1) is the radia- 
tion temperature of a blackbody at a temperature 
T, and Tb g is the cosmic background temperature 
of 2.7 K. The column density is then given by: 



Ntnt. = 



8ttv 3 Q(T ex ) e E -/ T ° 
c 3 g u A u i e h,/ / kT - 



1! 



rdv (5) 
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where Q(T ex ) is the partition function: 

oo 

Q(T ox ) = j^(2I + l)cxp(-Ei/kT cx ) (6) 

i=0 

In the case of the H 2 D + transition, g u = 9, A u i 
= 1.04 1CT 4 s"\ E u i = 17.9 K; in the case of the 
D 2 H+ transition, g„ = 9, k ui = 4.55 1CT 4 s _1 , E ut 
= 33.2 K. Using equation 4, we can estimate the 
upper limit on the D 2 H+ main beam temperature, 
under the assumption that the conditions of LTE 
are valid. With an excitation temperature of 8 
K, the maximum main beam temperature reaches 
0.53 K, which is about the observed main beam 
temperature of the D 2 H + tentative detection, in 
the case where the spatial extent is larger than the 
beamsize of 11". The derived column densities de- 
pend on the assumed value of the excitation tem- 
perature. Within the 7-9 K temperature range, 
the column density at the dust peak can vary by 
a factor of 2. At larger distance, considering the 
increasing kinetic temperature and the decreasing 
molecular hydrogen density, the excitation tem- 
perature could be as low as 6 K. However, at these 
positions, the H 2 D+ column densities should only 
be decreased by a factor ~ 2.5. Consequently we 
used a constant T ex for our observations as an 
approximation. The corresponding line parame- 
ters and column densities are presented in Table 1. 
The upper limit on the para-D 2 H + column density 
has been calculated using the thermal linewidth at 
8 K (0.27 km s" 1 ). 

At the (0",0") position, the column densities of 
ortho-H 2 D+ and para~D 2 H+ are 1.8 x 10 13 cm~ 2 
and <2.3 xlO 13 cm~ 2 , respectively (see Table 1). 
The abundances of species i, X{i), have been de- 
termined dividing the column densities N(i) by 
the associated H 2 column density derived from the 
1.3 mm dust continuum emission. At the dust 
peak, we obtain A(ortho - H 2 D+) ~ 1.5xl0" 10 
and A(para-D 2 H+) < 1.8xl0- 10 . 

We present in Figure 6 the variation of the 
observed ortho-H 2 D+, CO, H 2 , HCO+, DCO+, 
N 2 H + , and N 2 D+ column densities (crosses) 
across the core (as a function of the impact param- 
eter, the projected distance from the dust peak) 
as well as an upper limit on the para-D 2 H + . Note 
that we used the ortho-H 2 D + observations. We 
did not use any ortho-to-para ratio to estimate 
the total H 2 D+ column density or abundance be- 
cause of the large uncertainty on this ratio. A 



more thorough discussion on the orth-to-para 
ratio for H 2 D+ and D 2 H + will be performed in 
section 4.2.2. The observed column densities have 
then been averaged within the ranges delimited 
by the vertical dashed lines, at the positions (2i 
+ 1) x 10" where i=0,l,2,3... and are repre- 
sented by triangles. The same computation was 
performed to present in Figure 7 the variation of 
the abundances as a function of the distance to 
the core center, limited to the central 70" (r ~ 
9800 AU). We will compare in section 4.1 these 
observations with the result from a best-fit model 
(dashed lines). 

From Figures 6 and 7, we see that only CO is 
strongly depleted in the core center. Note that 
we used, in the CO column density computa- 
tion: 16 O/ 18 O=560 (Wilson & Rood 1994), and 
18 0/ 17 0=4 (Wouterloot et al. 2005; Ladd 2004). 
Defining the CO depletion factor, /d, as the ra- 
tio of the CO "canonical" abundance ([CO]/[H 2 ] 
= 9.5xl0~ 5 ; Frerking, Langer & Wilson 1982) 
and the observed CO abundance (N(CO)/N(H 2 )), 
Caselli et al. (1999) found /b = 10 toward the dust 
peak and concluded that the most likely explana- 
tion for the low CO abundance is the freeze-out 
of CO onto dust grains at high densities (n > 10 5 
cm -3 ). The corresponding radius of the depleted 
region is r - 6500 AU (~ 45"). 

The CO abundance is a critical parameter in 
the deutcration of the molecular ion H^" (see Fig- 
ure 1). In fact, from the abundance profiles pre- 
sented in Figure 7 it is clear that the degree of 
deuterium enhancement (with DCO + , N 2 D + , and 
H 2 D + ) increases toward the dust peak emission of 
L1544 where CO is highly depleted, as previously 
found by Caselli et al. (2002b). N 2 H+ does not 
show any signs of depletion. It is mainly formed 
by interaction between H 3 + and molecular nitro- 
gen, which is likely to be the main repository of 
nitrogen in the gas phase. N 2 is only slightly less 
volatile than CO (factor of ^0.9; Oberg et al. 
2005), so that the two neutrals are expected to 
behave similarly. However, N 2 H + is destroyed by 
CO (Bergin et al. 2001, 2002; Pagani et al. 2005; 
Aikawa et al. 2005), so that the CO freeze-out im- 
plies a drop in the destruction rate of N 2 H+ , which 
at least partially balance the lower formation rate 
due to the N 2 freeze-out (see also Aikawa et al. 
2001 for a discussion on this point). In fact, N 2 H+ 
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is observed to survive in the gas phase at higher 
densities (~10 6 cm~ 3 ) compared to CO (~ 10 5 
cm -3 ). This is also confirmed by the deuterium 
fractionation observed in N 2 H + (~0.2), about 5 
times larger than that measured in HCO + (Caselli 
ct al. 2002b). HCO+ is mainly formed via H| + 
CO and destroyed by dissociative recombination, 
so that its abundance is simply reduced by the 
freeze-out of CO, its parent species. From Figure 
7 it seems that the HCO + abundance is reduced 
at the dust peak, and increases at larger distance. 

3.1. Correlations 

In Figure 8, we show the correlation between 
the ortho-H 2 D + abundances at the 0", ± 20" and 
± 40" distance from the dust peak and the CO 
depletion factor, the DCO + /HCO + ratio and the 
N 2 D + /N 2 H + ratio. As intuitively expected, the 
ortho-H 2 D + abundance appears to be well corre- 
lated with the CO depletion factor (see Figure 1). 
Also, the degree of deuterium enhancement in the 
HCO + and N 2 H + molecules (measured from the 
DCO+ /HCO+ and N 2 D+ /N 2 H+ ratios) increases 
linearly with the ortho-H 2 D + abundance. The \ 2 
parameter is calculated for the three points where 
the impact parameters are 0", 20" and 40": 

x 2 = ^f WO-XflA 2 (7) 

where X f, s and X fa are the observed and fit values 
of the abundance respectively, and <Jx obs is the un- 
certainty in X b s . The associated probabilities are 
reported when a correlation is established. The 
surprisingly high confidence level for the correla- 
tions between H 2 D+ and the degree of deutera- 
tion in the HCO+ and N 2 H+ molecules (~ 100%), 
confirms that H 2 D + dominates the fractionation 
of these molecules at low temperatures. 

4. Chemical modelling 

In this section we will interpret the observations 
using chemical models. We will adopt a two-way 
strategy. First (in section 4.1), we use a full chem- 
ical model applied to a density structure derived 
from continuum observations to produce an over- 
all fit to all line observations presented in the pre- 
vious section. In a second step (section 4.2), we 



use a simplified chemical model focussing on the 
chemistry of H^" deuteration, in order to better 
understand the relation between CO depletion and 
deuteration, and even to provide some estimates 
of the ortho-to-para ratio in the deuterated forms 
of Hjj" that can be derived from our observations. 

4.1. The "best-fit" model 

To more quantitatively understand the column 
density and abundance correlations shown in Fig- 
ures 6, 7 we used the model described in Caselli 
ct al. (2002b), updated so that it now includes 
the multiply deuterated forms of H3 (as in Crapsi 
et al. 2005) and new values of the binding ener- 
gies of CO and N 2 , following the measurements 
by Oberg et al. (2005) as well as other modifi- 
cations to better account for the physical struc- 
ture of the core. Briefly, the model considers a 
spherically symmetric cloud, with the density gra- 
dient as derived by Tafalla et al. (2002), using 
the 1.3 mm dust continuum emission data from 
Ward-Thompson et al. (1999), where the cen- 
tral density is n(H 2 ) = 1.4xl0 6 cm~ 3 (see Figure 
5). The temperature profile has been included, 
using the recent findings of Young et al. (2004), 
where the temperature is about 7.5 K at the cen- 
ter and reaches about 12 K at the edge (see Fig- 
ure 5). The chemical network contains CO, O 
and N 2 , which can freeze-out onto dust grains 
and desorb due to cosmic-ray impulsive heating 
(as in Hasegawa and Herbst 1993). The initial 
abundances are: JQ(CO)=9.5xlO -5 (cf. Frerking 
ct al. 1982), A,(N 2 )=4.0xl0- 5 (slightly smaller 
than 6.6xl0 -5 , the cosmic value from Snow & 
Witt 1996, assuming that all nitrogen is in molec- 
ular form), and X i (0)=A i (CO)/2, a factor of two 
less than what is typically found in gas-phase-only 
chemical models (e.g. Lee et al. 1996). The abun- 
dances of the molecular ions (N 2 H + , Ii30 + and 
HCO + and their deuterated forms) are calculated 
in terms of the instantaneous abundances of neu- 
tral species, assuming that the timescale for ion 
chemistry is much shorter than that for freeze- 
out. The electron fraction has been computed, 
as described in Caselli et al. (2002b), using a 
simplified version of the reaction scheme of Ume- 
bayashi & Nakano (1990), where the chemistry of 
a generic molecular ion "mH+" is taken into ac- 
count, assuming formation due to proton transfer 
with H and destruction by dissociative recombi- 
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nation with electrons and recombination on grain 
surfaces (using rates from Draine & Sutin 1987). 
The "MRN" grain size distribution has been used 
(Mathis, Rumpl & Nordsieck 1977). We adopted 
the rate coefficients for the proton-deuteron ex- 
change reactions recently measured by Gerlich et 
al. (2002). The model is run until the column 
density of C 17 toward the core center reaches the 
observed value (N(C 17 0) = 6xl0 14 cuT 2 ; Caselli 
et al. 2002b). 

The "best-fit" parameters of the model, which 
best reproduce the observed molecular column 
densities at the dust peak and the observed col- 
umn density and abundance profiles, are the fol- 
lowing: 

• a cosmic-ray ionization rate of £ = 1 . 3 x 1 ~ 1 ' 
s _1 , standard value, typically used in chemi- 
cal models (since Herbst & Klemperer 1973), 

• binding energies for CO and N2 of -Ed (CO) 
= 1100 K and E D (N 2 ) = 900 K, values close 
to the binding energies measured for CO on 
water ice (Collings et al. 2003); the ratio 
between N 2 and CO binding energies (0.8) 
being comparable to the value (0.9) recom- 
mended by Oberg et al. (2005), 

• a binding energy for atomic oxygen of Ed (O) 
= 750 K, used in current gas-grain chemical 
models (e.g. Aikawa et al. 2005), 

• a minimum size of the dust grains, in the 
MRN distribution of a min = 5xl0~ 6 cm (10 
times larger than in MRN), 

• a sticking coefficient of S = 1.0 (Burke & 
Hollenbach 1983), 

• an initial abundance of metals (assumed to 
freeze-out with a rate similar to that of CO), 
X(M+) = 10~ 6 . 

In Figures 6, 7, 8 we overlaid the results from 
the best-fit model (in dashed lines), to the ob- 
servations. For the H 2 D + plots, we present the 
ortho-H 2 D+ observation points, and scaled the 
total-H 2 D + result from the model by 2.3, the fac- 
tor between the predicted and the observed value 
at the dust peak position (assumed constant across 
the core). From this, an ortho-to-para ratio of ~ 
0.8 can be deduced for H 2 D+, but, considering the 



uncertainties associated with this parameter (e.g. 
Pagani et al. 1992, Flower et al. 2005), we will 
postpone a discussion on this topic in the following 
section, where a parameter space exploration will 
be presented. For the D 2 H+ plot, we present the 
upper limit on the observed para-D 2 H + compared 
with the (total) D 2 H + result from the model. The 
ortho-H 2 D+ column densities and abundances ob- 
served at 0" , ± 20" and ± 40" are well reproduced 
by the model, within the error bars (see Figures 6, 
7). Note that although we need to assume a high 
degree of CO depletion in order to explain the CO 
observations, it appears that the HCO + column 
density is slightly under-estimated for the model 
at the dust peak emission. However, the strong 
variation (factor of 2) seen in Figure 6 between 
and 15" could decrease the central HCO + column 
density. The discrepancy between the HCO + and 
DCO + column densities at larger distances can 
be explained by the uncertainties on the opti- 
cal depth, since the less optically thick isotopes 
(HC 18 0+ and D 13 CO+) have only been used for 
the central position. Also the beamsize for the 
HC 18 + observations is larger (by 50%) than the 
20" range. 

Some N 2 depletion is needed in order to explain 
the N 2 H + and N 2 D + observations. Through the 
hyperfine structure of these species we can deter- 
mine directly the optical depth in several transi- 
tions using the relative intensities of the hyperfine 
satellites. This considerably reduces the error in 
our computations, compared to other species like 
HCO+ and DCO+. 

This detailed model of the ion chemistry in 
L1544 simulates the observed depletion in the core 
center, and can reproduce the observed depen- 
dency of the column densities and abundances of 
species such as N 2 H+, N 2 D+, HCO+, DCO+ as 
a function of the impact parameter. This allows 
us to separately discuss the relative contributions 
from the high-density depleted core and the lower 
density foreground (and background) gas. 

4.2. Chemical parameter space explo- 
ration 

In order to focus and analyse in detail the H 2 D + 
chemistry, we performed a parameter study using 
a model that computes the deuterated forms of 
H3 as a function of some key parameters, like 
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the grain size, the age of the L1544 condensation, 
and the cosmic rays ionization rate. This method 
has the advantage of concentrating on the H 2 D + 
chemistry, avoiding to reproduce other molecular 
observations. Before discussing the comparison of 
the theoretical predictions with the observations, 
we give a short description of the model used. 
It is an adaptation of the Ceccarelli & Dominik 
(2005) model, which has been developed for the 
proto-planetary disks. It computes the Hj, H 2 D + , 
D 2 H+ and Dg abundances in cold and dense gas. 
Since the involved temperatures (< 30 K) and 
densities (> 10 5 cm~ 3 ) are very similar to those 
found around L1544, the model can be used di- 
rectly, by just changing the geometry. For an easy 
and straightforward comparison with the obser- 
vations we compute the Hg chemistry in a gas 
cube with a given density and temperature. The 
relative abundances of the Hj deuterated forms 
are computed solving the charge equilibrium equa- 
tions and the deuterium chemistry equations. 
In this model we also consider grains as a possible 
source of Hg , H 2 D+, D 2 H+ and Dg destruction. 
In practice, the larger the CO depletion, the larger 
the H 2 D+ /H+ and D 2 H+ /H 2 D+ ratios. 
Two factors (other than the dust temperature) 
can modify the CO depletion : the age of the con- 
densation (larger ages give larger CO depictions 
because the molecules have more time to freeze out 
onto the grains), and the gas density (the conden- 
sation rate is proportional to the gas density). In 
addition, the cosmic ray ionization rate regulates 
the overall ionization degree in the condensation, 
and therefore the H;j" isotopomers abundances. 
Finally, the grain sizes enter both in the CO con- 
densation rate (via the grains area), and in the 
charge balance, because negatively charged grains 
recombine with the positively charged molecular 
ions. In this model, we adopted the same pa- 
rameters (binding energy for CO and N 2 , sticking 
coefficient) chosen for our best-fit model (see sec- 
tion 4.1). 

4-2.1. H 2 D + and D 2 H + versus CO depletion 

In Figure 9, we present the results of the model, 
varying the three key parameters of the model 
(cosmic ionization rate, the grain radius and the 
age of the core) in order to reproduce the total 
(ortho + para) H 2 D + and D 2 H + abundances. We 



plot the abundances as a function of depletion be- 
cause this parameter is more directly observed (via 
CO column density and dust continuum observa- 
tions) than the gas density. We also make a com- 
parison with the observations of ortho-H 2 D+ and 
para-D 2 H + to get an insight into the deuterium 
chemistry of Hg^ . 

We fixed the temperature of the cloud in the model 
to 8 K, which is within the range found from 
molecular and dust measurements. For compar- 
ison we also ran the cases with a temperature of 
10 K and did not find any substantial difference. 
The free-fall time is approximately 3 x 10 4 years 
for a density of about 10 6 cm~ 3 . In presence of a 
magnetic field, the collapse time is about an order 
of magnitude larger (Ciolek & Basu 2000). We will 
vary the age between 10 4 to 10 6 years, the latter 
being close to the lifetime of a starless core. In our 
calculations, we assumed all the grains in the cloud 
to have the same size, but we investigated the re- 
sult for different values of the grain radius. A grain 
size of 0.1 ^tm is the typical value assumed in chem- 
ical models for the interstellar medium, where it 
follows the MRN distribution. In the upper plot of 
Figure 9 we used a standard cosmic ionisation rate 
of 3 10~ 17 s _1 , a typical age of 10 5 years and var- 
ied the dust grain averaged sizes between 0.05 ^m 
and 0.2 ^m. In the middle plot we fixed the age to 
10 5 years, the grain size to 0.1 /zm and varied the 
cosmic ionization rate between 3 10~ 18 to 3 10~ 16 
s . In the lower plot we fixed the grain size to 0.1 
/xm, the cosmic ionization rate to 3 10~ 17 s _1 and 
varied the age of the cloud between 10 4 and 10 6 
years. The observation points (or upper limit in 
the case of D 2 H + ) and their corresponding error 
bars are supperposed in these plots: ortho-H 2 D+ 
on the left side and para-D 2 H + on the right side. 

As the grain size increases while the total grain 
mass is conserved, the abundance of grains relative 
to H 2 decreases oc a gr 3 in , and also the grain surface 
area per H 2 decreases oc GC ra ; n - As this effect slows 
down the freezeout of CO, the same CO depletion 
is therefore reached either after a longer time, or at 
a higer density. Since we keep the age constant, 
the effect of the density is observed in Figure 9: 
larger grain sizes correspond to higher densities, 
at which the overall degree of ionization is smaller. 
Since H 2 D + is the dominating ion, this is directly 
mirrored in the H 2 D + abundance. Also, decreas- 
ing the cosmic ionization rate leads to a decrease 
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in the abundances. Indeed, H3" 1 " ions (and conse- 
quently their deuterated forms) are formed by the 
ionization of H2 due to cosmic rays. And finally, 
increasing the age of the cloud will increase their 
abundances, as the CO depletion rate is time de- 
pendant. Consequently, for a more evolved cloud 
the same CO depletion is achieved for lower den- 
sities, corresponding to a higher degree of ioniza- 
tion, which is again directly reflected in the H 2 D + 
and D 2 H + abundances. 

4-2.2. The ortho and para forms 

Both H 2 D + and D 2 H + molecules have ortho 
and para forms, corresponding to the spin states of 
the protons (for H 2 D + ) or deuterons (for D 2 H + ). 
In order to compare the modeled abundances with 
the observations of one spin state only, it is criti- 
cal to know the ortho-to-para ratio for these two 
molecules. Under LTE conditions, at temperature 
T, the relative populations of the lowest ortho 
(li,i) and para (0q,o) levels of H 2 D + would be: 



n(li.i) . , 86.4 
- 9 x exp( — ) 



(8) 



n(0 0l0 ) ' T 
and the relative populations of the lowest ortho 
(0o,o) and para (lo,i) levels of D 2 H+ would be: 



n(lo,i) 



9 . 50. 2 N 

- x expl 

6 T 



(9) 



n(0 ,o 

With these formulae, at 8 K, the H 2 D+ ortho- 
to-para ratio would be ~ 1.8 1CP 4 and the D 2 H+ 
para-to-ortho ratio would be ~ 2.8 10~ 3 . The 
ortho form of H 2 D + is produced mainly in reac- 
tions of the para form with ortho-H 2 (e.g. Gerlich, 
Herbst & Roueff 2002). Therefore, its high o/p ra- 
tio is attributable to the relatively high ortho-H 2 
abundance as first noticed by Pagani et al. (1992) . 
Because the o/p ratio is not thermalized at the low 
temperature considered here, the o/p H 2 D + ratio 
is not thermalized either. This can be illustrated 
in Flower, Pineau des Forets, Walmsley (2004) 
model where, at temperatures lower than 10 K, a 
hydrogen density of 2 x 10 6 cm~ 3 and a grain size 
of 0.1 /um, the o/p-H 2 D + reaches unity and the 
p/o-D 2 H + value is about 0.1. Increasing the grain 
size will decrease the grain surface, leading to a 
decrease of the H 2 formation rate. Therefore, the 
H 2 ortho-to-para ratio will obsiously decrease, as 
well as the H 2 D+ ortho-to-para ratio. 



From our observations we find that para- 
D 2 H+/ortho- H 2 D+ is less than 1.3 at the dust 
peak emission assumed to be at 8 K. In the prestel- 
lar core 16293E (Vastel et al. 2004) we measured 
a para-D 2 H + /ortho-H 2 D+ of 0.75 for an excita- 
tion temperature of 10 K. We also computed the 
H 2 D + ortho to-para ratio and an upper limit on 
the D 2 H + para-to-ortho ratio by comparing the 
observed abundances of ortho-H 2 D + and para- 
D 2 H + with the total (ortho + para) abundances 
of H 2 D+ as calculated using the model described 
in the previous section. In Table 2, the ortho-to- 
para ratio for H 2 D + and the para-to-ortho ratio 
for D 2 H + are quoted in order to reproduce the 
values obtained by the model for different sets of 
parameters, which are the cosmic ionization rate, 
the age of the core and the grain radius. We 
can directly compare our results on H 2 D + with 
the Flower, Pineau des Forets, Walmsley (2004) 
model, even if their study assumes complete deple- 
tion (i.e. that CO abundance should be less than 
10~ 6 ). Indeed the abundance of both ortho and 
para spin states of H 2 D + depends on the ortho 
and para forms of molecular hydrogen (through 
the proton-exchanging reaction of H 3 + with H 2 
followed by reaction 1) which does not vary as a 
function of depletion. On the contrary, the abun- 
dance of both para and ortho forms of D 2 H + is 
determined by reactions with HD (produced on 
the grain surfaces) and will therefore depend on 
the core depletion. The Flower et al. model pre- 
dicts H 2 D+ ortho-to-para ratios larger than the 
maximum value of 0.3 we observed for a temper- 
ature of 8 K and a H 2 density of 2 x 10 6 cm -3 
(Pineau des Forets, private communication) span- 
ning ranges up to 0.4 /im of the grain radius. As a 
consequence, since the H 2 D + ortho-to-para ratio 
decreases as a function of the grain radius, it is 
likely that this should be larger than 0.3 /im. This 
depletion of small grains in this core is consistent 
with grain coagulation since ice condensation is 
not enough to increase the grain radius. 

5. Conclusions and Perspectives 

In this paper we studied the prestellar core 
L1544, focusing on the H 2 D + chemistry through- 
out the cloud. It is now widely accepted that the 
H 2 D + molecule is the main tracer of the CO de- 
pleted prestellar cores, and we point out in this 
paper that the H 2 D + emission is extended (over 
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60"), and an excellent tracer of the dust contin- 
uum, with a emitting radius of about 7,000 AU 
in the case of LI 544. Hence, the line profile of 
this molecule would provide a crucial guide to the 
dynamical behavior of the high-density core. It 
is likely that the double-peak profile found in the 
central position, as well as positions around, is 
broadened by the central infall and is absorbed 
in the outer parts of the core (van der Tak et al. 
2005). 

We first used a model of the ion chemistry coupled 
with the physical structure of the core of L1544, 
including the deuterated isotopologues of the H3 4 " 
ion. This simulates the observed depletion and 
can approximatively reproduce the observed de- 
pendence of the column densities of species like 
N 2 H+, N 2 D+, HCO+, DCO+, H 2 D+ as a func- 
tion of radius. 

This study reveals a correlation between the 
ortho-H 2 D + abundance and 1) the CO deple- 
tion factor, 2) the degree of deuteration in the 
HCO + molecule, 3) the degree of deuteration in 
the N 2 H+ molecule. H 2 D + will survive longer, at 
higher density than N 2 H+ and N 2 D+. 
We then used a simpler model focusing on the 
H 2 D + chemistry where we did a wide parameter 
study. We discuss how the H 2 D + and D 2 H+ abun- 
dances depend on the cosmic ionization rate, the 
age of the core, and the grain radius, by varying 
these parameters. It appears that the most impor- 
tant parameters to reproduce the observed values 
is the grain radius as small grains accelerate the 
freeze-out of CO and the observed values are con- 
sistent with a freeze-out rate dominated by larger 
grains. Therefore, we found that to reproduce the 
observations we need a larger grain radius of 0.3 
fim. 

This study can be considered as a springboard for 
observations to come, since the current observa- 
tions are limited by the poor atmospheric trans- 
mission at the relevant frequencies. Table 3 lists 
some of the major telescopes and interferometers 
that can be used for the study of H 2 D + chemistry 
in prestellar cores, proto-planetary disks and pro- 
tostars. Probably, D 3 + can not be observable be- 
cause enhanced D3 + abundance implies very cold 
and very dense regions. Since D3 + is a symmetric 
molecule, it does not have rotational transitions 
and does have its bending modes are in the near 
infrared. Therefore, these transitions will only 



be observable in absorption against a strong near 
infrared continuum. H 2 D + and D 2 H + are hence 
the only tracers of cold, dense phases of molecular 
clouds prior to star formation. 
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Fig. 1 . — Main reactions involving the deuterated 
forms of the molecule. When CO and N2 are 
depleted, the molecular reactions presented with 
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Fig. 2. — Ortho-H2D + and para-D2H + observa- 
tions at the dust peak emission corresponding to 
the (0",0") position. The temperatures are T* 
in Kelvins. The vertical solid line corresponds to 
the velocity center of the D 2 H + gaussian fit (7.29 
km s _1 ). 
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Fig. 3. — Map of the H 2 D+ (lio-ln) line centered 
on the dust peak of L1544. The Y-axis represents 
the main beam temperature. The position is indi- 
cated in arcseconds in right ascension and declina- 
tion offsets from the central position, on the top 
left of each spectrum. The reference for the veloc- 
ity center at the (0",0") position (7.28 km s _1 ) is 
indicated in dashed lines. 
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Fig. 4. — Integrated intensity maps of H 2 D+ (lio- 
lu), N 2 H+ (1-0) and N 2 D+ (2-1) superposed on 
the 1.3 mm continuum emission map smoothed to 
a resolution of 22" (gray scale). Contour levels 
are 30%, 50%, 70% and 90% of the peak (0.54 
K km s- 1 for H 2 D+), and 50% of the peak (5.5 
K km s- 1 for N 2 H+ and 2.1 K km s" 1 for N 2 D+). 
The observed positions in H 2 D + are reported as 
triangles. The (0",0") position corresponds to 
a 2 ooo = 05 71 04 m 17 s .21, <5 2000 = +25° 10' 42.8". 
The beam sizes of the observations are also re- 
ported in the lower right corner of the figure. 
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Fig. 5. — Density (plain line) and temperature 
(dashed line) profiles as a function of the radius 
(in astronomical units) or the impact parameter 
(in arcseconds), used for the best-fit model. 
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Fig. 6. — Variation of the CO (in a 20" beam), 
H 2 , HCO+, DCO+, N 2 H+, N 2 D+, H 2 D+ and 
D 2 H+ column densities as a function of distance. 
The cross represent the observation points and 
the 3 a error, the triangles represent the observa- 
tion points averaged in the bins delimited by the 
dashed lines, and the dot-dashed lines represent 
the results from a "best-fit" model. For H 2 D+, 
note that the variation of the observed ortho col- 
umn density is compared with the modeled ortho 
variation (see text). Also for D 2 H+, the upper 
limit on the para column density is reported on 
the plot of the modeled ortho + para variation. 
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Fig. 7. — Variation of the CO (in a 20" beam), 
HCO+, DCO+, N 2 H+, N 2 D+, H 2 D+ and D 2 H+ 
abundances as a function of distance. The tri- 
angles represent the observation points averaged 
in the bins delimited by the dashed lines, and 
the dot-dashed lines represent the results from a 
"best-fit" model. For H 2 D + , note that the varia- 
tion of the observed ortho abundance is compared 
with the modeled ortho variation (see text). Also 
for D 2 H + , the upper limit on the para abundance 
is reported on the plot of the modeled ortho + 
para variation. 
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Fig. 8. — Variation of the observed ortho-H 2 D + as 
a function of the depletion factor, DCO + /HCO + 
ratio and N 2 D + /N 2 H + ratio. The x 2 parameter 
and its probability are reported when a correlation 
is established. Superposed are the results from the 
best-fit model (see section 4.1) 
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Fig. 9. — Variation of the H 2 D + abundance (re- 
spectively D 2 H+) for the model (plain line) as 
a function of CO depletion factor. The points 
with the corresponding error bars represent the 
observed abundances of ortho-H 2 D+ and para- 
D 2 H+ towards LI 544 whereas the lines show the 
modeled abundances of the ortho + para states. 
Hence we expect the observed values to lie below 
the modeled ones. The temperature was fixed to 8 
K. In the upper plot we varied the grain size: 0.05, 
0.1, 0.2, 0.3 and 0.4 /im. Note that increasing the 
grain size decreases the H 2 D+ and D 2 H + abun- 
dances (see text). In the middle plot we varies 
the cosmic ionization ray: 3 10" 16 (plain line), 3 
10" 17 (dashed line) and 3 10" 18 s" 1 (dot-dashed 
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Table 1: Lines parameters, and column densities 
of the ortho-H2D + and para-D2H + observations 
for an excitation temperature of 8 K. The 1 a er- 
rors arc noted in parenthesis. The upper limit (3 
a) on the para-D 2 H + column density has been 
determined using a linewidth of 0.27 km s _1 (see 
section 3). 



Position 


rms 




At; rcs 


Av 


J'T mh Av 


N(o-H 2 D+) 


tau 


(") 


(K) 


km s _1 


km s _1 


km s -1 


mK km s _1 


cm -2 




(0,0) 


0.11 


7.277(0.016) 


0.039 


0.50(0.03) 


495(33) 


1.78(0.10) 10 la 


0.58(0.05) 


(20,0) 


0.12 


7.319(0.045) 


0.039 


0.52(0.12) 


255(37) 


7.77(1.01) 10 12 


0.24(0.04) 


(0,-20) 


0.10 


7.382(0.040) 


0.039 


0.59(0.09) 


232(33) 


6.92(0.90) 10 12 


0.19(0.03) 


(-20,0) 


0.13 


7.269(0.033) 


0.039 


0.52(0.06) 


293(40) 


9.14(1.10) 10 12 


0.29(0.05) 


(-20,-20) 


0.11 


7.238(0.023) 


0.039 


0.20(0.06) 


98(20 ) 


2.99(0.58) 10 12 


0.24(0.05) 


(-40,0) 


0.13 




0.077 




< 49 


< 1.46 10 12 




(-40,20) 


0.12 


7.294(0.035) 


0.039 


0.27(0.07) 


115(26) 


3.45(0.73) 10 12 


0.21(0.05) 


(-20,20) 


0.21 


7.316(0.018) 


0.039 


0.37(0.04) 


386(54) 


1.28(0.15) 10 13 


0.58(0.10) 


(0,20) 


0.08 


7.094(0.015) 


0.039 


0.43(0.04) 


288(22) 


9.27(0.62) 10 12 


0.35(0.03) 


(20,20) 


0.09 


7.256(0.026) 


0.039 


0.29(0.08) 


111(15) 


3.30(0.57) 10 12 


0.19(0.03) 


(0,40) 


0.11 


7.300(0.024) 


0.039 


0.18(0.05) 


75(20) 


2.24(0.55) 10 12 


0.20(0.03) 


(-20,40) 


0.12 


7.138(0.032) 


0.039 


0.45(0.09) 


206(34) 


6.24(0.94) 10 12 


0.23(0.04) 


(-40,40) 


0.10 




0.039 




< 27 


< 4.52 10 12 




(20,-20) 


0.14 


7.378(0.040) 


0.039 


0.49(0.06) 


224(42) 


6.79(1.15) 10 12 


0.23(0.04) 


(20,40) 


0.12 




0.077 




< 46 


< 1.27 10 12 




(20,-40) 


0.11 


7.436(0.054) 


0.077 


0.46(0.10) 


215(45) 


6.55(1.23) 10 12 


0.23(0.04) 


(40,-40) 


0.12 




0.077 




<46 


< 1.27 10 12 




Position 


rms 




Av ICS 


Av 


/T mh Av 


N(p-D 2 H+) 


tau 


(0,0) 


0.072 


7.290(0.007) 


0.021 


0.27 


< 103 


< 2.27 10 la 


< 1.13 
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Table 2 

h 2 d+ ortho-to-para ratio and upper limit on the d 2 h+ para-to-ortho ratio, varying the 
cosmic ionization rate, the core age and the grain radius. 





parameters 




- H 2 D+ 


E D 2 H+ 

o A 


cosmic ionization rate age 
(s^ 1 ) (years) 


grain radius 

(H 


3 10- 17 


10 5 


0.05 


0.02 


< 1.07 


3 10- 17 


10 5 


0.1 


0.03 


< 0.15 


3 10- 17 


10 5 


0.2 


0.08 


< 0.55 


3 lcr 16 


10 5 


0.1 


0.03 


< 0.29 


3 10- 18 


10 5 


0.1 


0.07 


< 0.29 


3 10- 17 


10 4 


0.1 


0.07 


< 0.29 


3 10- 17 


10 6 


0.1 


0.03 


< 0.29 


3 10- 17 


10 5 


0.3 


0.17 


< 0.63 


3 10- 17 


10 5 


0.4 


0.32 


No solution 


3 10- 18 


10 5 


0.4 


0.33 


< 11.81 


3 10~ 17 


10 4 


0.4 


0.33 


< 11.81 
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Table 3 

Current and future facilities for the chemistry of H 2 D+ and D 2 H+ in prestellar cores, 

PROTO-PLANETARY DISKS AND PROTOSTARS. 



Name 


Aperture 


Platform 


Available 




H 2 D+ 




D 2 H+ 










ll,0 _ ll,l 


lo,i-0o,o 


li,o _ lo,i 


li,i-0o,o 










(372.4 GHz) 


(1.37 THz) 


(691.7 GHz) 


(1.48 THz) 


CSO 


10.4 m 


Mauna Kea (USA) 


Y 


Y 


N 


Y 


N 


JCMT 


15 m 


Mauna Kea (USA) 


Y 


HARP B 


N 


Y 


N 


SOFIA 


2.5 m 


Airborne (747) 


2007 


N 


Casimir 


Casimir 


Casimir 










GREAT (CONDOR) 




GREAT (CONDOR) 


Herschcl (HIFI) 


3.5 m 


Space (L2) 


2007 


N 


N 


Y 


Y 


ALMA 


50 X 12 m 


Atacama (Chile) 


2010 


Y 


N 


Y 


N 


APEX 


12 m 


Atacama (Chile) 


2005-2006 


Y 


CONDOR 


Y 


CONDOR 



